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Abstract 

For the first time, we obtain the entropy variation law 
in a laser process after finding the Kraus operator of 
the master equation describing the laser process with 
the use of the entangled state representation. The 
behavior of entropy is determined by the competition 
of the gain and damping in the laser process. The 
photon number evolution formula is also obtained. 



1 Introduction 

Since the theoretical foundation proposed by Albert 
Einstein in 1917 and the building of first function- 
ing laser by Theodore H. Maiman in 1960, laser has 
been successfully applied in various of areas, includ- 
ing laser cooling technique developed by Steven Chu 
et al [21 S. As one of the most important concept in 
physics, entropy measures the disorder of a system. 
Studying the evolution of the entropy, we can get 
a clear understanding of how a laser beam is created 
by appropriate pumping. A few works had been done 
concerning the entropy exchange between a laser and 
its environment 5113. However, the evolution of en- 
tropy in a laser itself has not yet been studied before. 
In this work we shall derive the entropy evolution 
law of a laser process. Our results explain how the 
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self-organization phenomenon happened in a laser. 

In quantum optics theory the time evolution of 
laser in the lowest-order approximation can be de- 
scribed bj^t_he_following master equation of density 
operator 
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where g and k are the cavity gain and the loss, re- 
spectively, , a are photon creation and annihilation 
operator, respectively. It is also known that the evo- 
lution due to the interaction between a system and 
its environment can be ascribed to an evolution from 
the initial density operator po to p (t) 
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such an expression is named an operator-sum (Kraus) 
representation, Af„ is named Kraus operator. So far 
as our knowledge is concerned, the entropy variation 
in laser-channel has not ever been reported. In this 
paper we shall show how the entropy of an initial co- 
herent state Po = \z) {z\ (the fact that n- photon dis- 
tribution in a coherent state is Poisson distribution 
exactly fits the measurement result of photon distri- 
bution in a laser light) varies in the laser process, 
before doing this, we should first derive the Kraus 
operator by solving the master equation ([l}. 

Our way is introducing the two-mode entangled 
state 
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where is a fictitious mode independent of the real where the n-mode coherent state is defined as 
mode , |0) is annihilated by a, [a,a^] — 
state \r] = 0) possesses the properties 
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Operating the both sides of ([1]) on the state [77 = 0) = 
|/) , and denoting \p) = p\I) , and using ^ we have 
the time-evolution equation for \p (t)) , 
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where |po) = Po |^) i Po is the initial density operator. 
The formal solution of ([5]) is 
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and 



U (t) = exp 
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It challenges us how to disentangle the exponential 
operator U (t) . This reminds us of two theorems 
about the normally ordered expansion of multimode 
bosonic exponential operators, which is helpful to dis- 
entangle U (t). 

2 Two Theorems 

In order to find the disentangled form of ([7]) we 
employ two new theorems about the normally or- 
dered expansion of multimode bosonic exponential 
operators : 

Theorem 1: The multimode bosonic exponential 
operator expH, where H =^BTB, B is defined by 
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r is a 2?i X 2n matrix, has its n-mode coherent state 
representation: 
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/„ is the n X n unit matrix. Q, L, N, P are n x n 
complex matrices 



matrix, obeying 
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or 



QL = LQ, QP-LN = I, (13) 
NP = PN, PQ~NL = I. 

Theorem 2: By performing the integration in © 
with the technique of integi'ation within an ordered 
product of operators , we have 

e^pn = ^eM-^,A\LP-^)An 
xexp{A'<{\nP-^)A}exp{^A{p-^N)A}. 



Now we first appeal to Theorem 1, so we should 
identify U {t) in ((T)) as expH, where A = ( a a ). 
After putting U (<) into the following symmetrized 
matrix form 
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with r being the symmetric matrix 
r = t 

here 
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we then follow PT|) to calculate exp(rn) with 

rn = t 

therefore 
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Thus according to Theorems 1 and 2 we have 
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where we have used 
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and using (|4]), ([6]) becomes 
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where 
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is the Kraus operator, and one can check 
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The Photon Number Evolu- 
tion 



Now we have an explicit solution of the density ma- 
trix of a laser (I^Bl) . we first calculate the evolu- 
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tion of photon number in a laser process initially 
in a pure coherent state, i.e., po = k) (-^l, \z) = 
exp[-|zp/2 -I- 7ntl 
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Then using |0) (0| —: e :, the normal or- 

dering form of the vacuum projector, as well as 
/ ^ I2) (zj = 1 we have the expected photon number 
evolution formula 
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(25) We can easily write down the asymptotic behavior 
of (n) when t — ^ -|-cx3 as the following: 

(26) If — 5; then (n) — \z\^ + '2gt , the photon number 
increases linearly with time. 
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li K < g, then (n) - + \z\^j e^^s-'^)t^ the 

photon number increases exponentially when t 
+00. 

the expected photon 



If K > g, then (n) 



number approaches a constant when t 
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4 The Entropy Evolution in a 
Laser 

We now calculate how the entropy of a laser evolves 
with time. Using (|26p , the density matrix p (t) of a 
laser initially in a pure coherent state \z) {z\ is 

p (t) = T3 exp e2(9-«)* IngTi 
2 : atiaJe^°*+^*'^~'^'^'^ : e'^^'^i"^^ (31) 
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By the Baker-Campbell-HausdorfF formula, if 
[X,Y]=\Y + ^l, 

then 
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we can compact the three exponentials in pip into a 
single exponential 
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with direct calculations. Thus we see how a pure 
state \z) {z\ evolves into a mixed state, so the entan- 
gled state representation in ^ can well expose the 
entanglement between the system and its environ- 
ment. Then the logarithm of p {t) can be evaluated 
as 
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Therefore the von- Neumann entropy of p {t) is 
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where fc^ is the Boltzmann constant. Since 
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Finally the entropy variation law is 
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Figure 1: S (p {t)) /fcs for ^ = 4, k = 1 and 5 = 2, 1, 0.5 
respectively 
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Figure 2: S{p{t)) /{{n)kB) for 2 = 4, k = 1 and g = 
2, 1, 0.5 respectively 

We also write down the asymptotic behavior of 
S {p (t)) when t — > +00 as the following: 

U K = g, then S (p (t)) /fc^ ~ 1 + In {2gt) as t 
+00, the entropy increases logarithmically. 

If K < 5 then S'(p(i)) /ks l + ln^ + 2 (g - n)t 
as t — > +00, the entropy increases linearly. 

If K > g then S (p {t)) /fc^ ^ In ^ + In | as 
t — >■ +00, the entropy approaches a constant. 

The results of expected photon number and en- 
tropy of the laser do not depend on the phase of pa- 
rameter z, as one should expect, since the absolute 
phase of z in a coherent state is non-physical. It is re- 
markable that the entropy is completely independent 
of z. 

Plots of S {p {t)) and the "specific entropy" ^^f^ 
in unit of ks for z = 4 (16 photons in average), k — 1 
and g = 2, 1,0.5 respectively are shown in Figure 1 
and 2. Besides the We can see clearly from the two 
figures that when the pumping rate g is less than the 
loss rate k, the photon number and entropy will ap- 
proach to constants, the photons are in fact in sort of 
thermo-equilibrium with an equivalent temperature 
T = 1^ In ^. When g is larger than k, while the en- 
tropy increases linearly with time, the expected num- 
ber of photons increases much more fast, therefore the 
specific entropy will goes to zero exponentially. The 
photons in the laser are highly coherent in this case. 
The above results indeicatc that a laser can gener- 
ate laser beam only if it works with sufficiently high 
pumping rate g. 
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